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^ ■ Abstract 

o 

■ The spin wave modes of cylindrical nanowires of moderate diameter-to-length ratio are investigated 

in this article. Based on three dimensional simulations and analytical calculations we determine the 
spatial structure of the modes. We show that standing spin waves and localized edge modes form the 
discrete spectrum of the nanowires. Using a simple analytical model we infer an extended dispersion 
relation for spin waves in cylinders. Considering the variation of the demagnetizing (internal) field we 
O . show that the localized dipole-exchange modes at the edges are always present. 

s 

^ Magnetic nanostructures are intensively studied nowadays for fundamental and practical purposes [T]. 
Understanding the magnetism in low-dimensions presents a fundamental interest, while achieving smaller 
I \ data storage media and magnetic memory is of utmost technological importance. Large periodic arrays of 
^ ■ magnetic nanowires have been easily produced with inexpensive techniques like electrodeposition[2j. For 
^ such structures, the manipulation of magnetization in very short times is very important for high-speed 
Q ^ applications and requires a complete understanding of the nature of magnetic excitations, magnons or 

spin-waves, and their dependence on geometry. 
^ . The spin- wave spectrum for ellipsoidal samples is known for some time[3[|l]. In uniform magnetized 

^ ' ellipsoidal samples the demagnetizing field is uniform. As many of the magnetic nanostructures studied 
nowadays have in general a non ellipsoidal form, the demagnetizing field inside them is nonuniform even 
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1 Introduction 



OO ■ if they are uniformly magnetized [5j. Thus the shape can drastically affect the dynamic properties and the 
^ ' spectrum of the spin- waves will be modified (see Ref.[6j and references therein). As the dimension of the 
wires decreases, the exchange interaction becomes important. We consider here that both dipolar and 
exchange interactions contribute to the spin-wave spectrum. Obtaining a general theory is challenging, 
but approximate analytical solutions and numerical results can be derived for certain geometries. 

Reducing the dimensions of magnetic nanostructures is the natural trend to obtain smaller devices. 
Until now, only very long nanowires were studied with diameter-to-length ratio (d/L) <C 1. In this paper, 
we investigate nanowires of moderate aspect ratio magnetized along the axial direction. We develop the 
pi I analytical theory and present numerical results of the spatial distribution of spin waves modes in this type 
of nanostructures. We will outline the important features of taking into account the in-homogeneity of the 
demagnetizing field for the confined geometry under study. The previous studies of cylindrical geometry 
neglected the nonuniformity of the demagnetizing fields considering the cylinders semi-infinite |i 



2 Analytical model 

We start with a ferromagnetic cylindrical sample, magnetized to saturation along the axial direction {z 
direction) which ensures that the longitudinal component of M is much larger than the transverses ones, 
M = {rrix^my, M^). A static magnetic field H is applied along z and an rf field h(r,t)=/i(r)e~^^^ is 
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applied in perpendicular direction, Hext — {hx^hy, H). We assume that the rf field and the saturation 
magnetization are uniform in the sample, even if the state of uniform magnetization is not actually 
realized in this geometry. Two main approaches are used for solving the spectrum of spin- waves in confined 
structures. In both, one solves simultaneously the linearized Landau-Lifschitz equation of motion for the 
magnetization together with the Maxwell equations satisfying the electromagnetic boundary conditions 
and the exchange boundary conditions [9J. The two methods are equivalent. We choose here the method 
of magnetic potential where the relation m(h) is found first from the equation of motion and then we 
search for the solutions of Maxwell equations which satisfy the boundary conditions. 

The Landau-Lifschitz (LL) equation of motion for the magnetization neglecting damping is: 

^(r,t) = -7MxHe// (1) 

where 7 is the gyromagnetic ratio. The effective field inside the sample represents the sum of the external 
applied field Hg^t, the exchange field iiexch — DV^Nl and the demagnetizing field Udemag — — NM 
(excluding crystal anisotropy): 

HgJJ = ^ext ~l~ ^exch ~^ ^demag (2) 

where D is the exchange stiffness and N the demagnetizing tensor field. For non-ellipsoidal bodies the 
demagnetizing tensor field (in the first order) is function of position and is defined in the Fourier space 
as in Ref.[10j: 

iV,,(k) = ^k,kj (3) 

with D(k) the shape function. The LL equation (Eq.([T])) describe an uniform precession of the magneti- 
zation around the effective field Heff- To find the normal modes of the magnetization we search for non 
zero solutions of the LL equation. Considering the dynamical magnetization uniform in the sample we 
find: 



iurrix — ^my[Hi — DV'^] — ^Mghy 
-iurriy ^mx[Hi - DV'^] - -fMshx (4) 

with Nzz and Nj^^ the longitudinal (dc) and transversal (ac) demagnetizing factors and H^=H-47r(N22- 
N^^)M5. Here, we'll consider as first order approximation that the demagnetizing tensor field can be 
diagonalized as A^22+2A^^^=47r. 

The demagnetizing field is determined by Maxwell equations in the magnetostatic limit: 



V X H,nt = (5) 
VB = V(H,,, + 4^M) = (6) 

where the internal field is Hint — h — Nrrm^ H — N^zMg. From Eq.([5j), the magnetic field may be written 
as h = — V^^rn, where is the scalar magnetic potential. Replacing the scalar potential in Eq.® we 
obtain: 

-V^«,„ + (4,-^„)(^ + ^)-M,^=0 (7) 

Eq.([7j) has the form of an anisotropic Laplace equation or Walker type equation for cylinders with the 
last term coming from the variation of the demagnetization (internal) field near edges. In the center of the 
cylinder, where the demagnetization field is almost constant (here the longitudinal demagnetization field 
is zero) the Walker equation reduces to the classical equation for cylinders [8j. We considered here only 
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Figure 1: Variation of the demagnetizing tensor components (here normahzed to Nzz^^^^rr—^) along the 
axial direction for a cylinder with length-diameter (L/d) ratio of 10: N^^ solid line, N^^ dashed line and 
dNzz/dz triangles. Only half of the cylinder is shown with largest variation near the edge. 

the variation of the longitudinal demagnetizing factor N^z along z (see Appendix for more details). The 
nonuniformity of Nzz strongly affects the mode frequencies as each mode feels a different demagnetizing 
factor. The full analysis should also take into account the variation of N^^ in the radial direction r, which 
we neglect here as Nrr is almost constant with r due to the fact that our cylinders have a reduced aspect 
ratio (diameter-to-length inferior of 0.1) and such an analysis will introduce a high degree of analytical 
complexity. 

We state that we made here the simplifying assumptions: the static magnetization is considered as 
constant and independent of the external field and that the dc and ac demagnetizing fields are the same 
and can be calculated with the expressions given in Eq.(|3j). Even if the non-diagonal elements of the 
demagnetizing tensor field are not zero, they can be overlooked as they are negligible (the N^^ component 
is inferior of 10~^). The variation of the demagnetizing factors [lOj along the axial direction is shown in 
the FiglH 

Choosing the scalar potential of the form ^^rn(^5 ^i z)=Jri(qr)exp(in0+ik2z) where Jn(^r) is the Bessel 
function of order n, we obtain the transcendental equation (see Appendix): 



X [2Hi + (47r - Nrr)Ms] + j\kl + q^)Hi [Hi + {4n - iV,,)M,] 

- fc2^2jy^,(47r - Nrr)Hi - J^MsD{4tT - Nrr)kl{kl + q^) 

- -^^^MsHf + -^2j^M,DHi - j^D^Ms^^ (8) 

In the limit of very long wavelength, when the wave vector has all components zero we derive the 
following expression: 



which takes the form of the Kittel uniform modeflj if we neglect the higher exchange terms (derivative 
of order three and five) and keep only the first derivative term: cjq — 7 Hi. However, no real uniform 
mode exists in cylinders of finite aspect ratio as the internal field varies near edges. In the case of a 
nonuniform magnetic field we usually assume that a spin wave can propagate with continuously changing 
wave vector kz{z)[12\. At certain points (surface), the wave vector becomes zero and the condition of 
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quasi-static variation of the internal field is no longer satisfied thus the solution of the equation of motion 
has to be found with variable parameters. This surface is called turning surface by analogy with quantum 
mechanics. The Eq.Q at the turning surface has no physical meaning. A localized spin wave is excited 
with changing wave vector propagating from the turning surface in the direction of the decreasing internal 
field. This effect called spin- wave well was demonstrated in stripes recently [13j. 

For samples of finite size, the dispersion relation is shape dependent. The complete spin- wave spectrum 
and the modes shape depend strongly on the boundary conditions. We expect some degree of pinning at 
the edges, and this degree of pinning is determined by a competition between dipolar magnetostatic energy 
and exchange energy. Apart from the Maxwell boundary conditions, continuity of the magnetization and 
the magnetic potential, usually a Rado-Weertmann type of boundary condition [HI [15] is considered at 
the cylinder surface: 



at z=±L/2 with L the length of the cylinder and p the so-called pinning parameter which is determined 
by the effective surface anisotropy K5 and the exchange stiffness constant D: p^K^/D. This condition 
implies stationary waves in the z direction with the particular solution of the type m^^Asin fc^z+Bcos k^z. 

The frequency of an eigenmode is constant throughout the magnet and the spin wave vector varies 
to accommodate the changes in the internal field profile. The dynamic magnetization has a plane- 
wave character. In almost all the volume of the cylinder, between the turning surfaces at both ends 
(excluding the edge domains) , the variation of the demagnetization field is small and usually the averaged 
value of the components of demagnetizing field is used. The derivatives of the N^^ in the Eq.(|8]) are 
negligible. For moderate aspect ratio cylinders we can use a simple and intuitive model where the 
dynamical magnetization takes the form: 



where is the longitudinal wave number for the Ith longitudinal mode (backward geometry) and n 
indexes the radial modes. This relation is valid between the turning surfaces, where the internal field is 
almost constant. In this region, stationary waves on the axial and radial directions are formed and we 
can have mixed modes. The distance between the turning surfaces, or effective length Az, is considered 
to be approximately 0.8L (with L the length of the cylinder). The dispersion curves are characterized by 
two quantization numbers n and /. The quantization parameter / should be taken with care, it provides 
here a qualitative description of longitudinal modes which is to be compared with the micromagnetic 
simulations. We can use the mean value of the wave vector of each mode which can be evaluated as: 
k^ = /tt/Az, where Az is the effective length where the mode is localized in the sample. The spectrum 
of the SW consists of a series of dispersion curves characterized by two quantization numbers n and I 
and depend on the ratio of the cylinder dimensions. Changing the aspect ratio of the cylinders induces a 
change in the values of wavenumbers (and also a change in the demagnetizing field). For example, if the 
aspect ratio is diminished, q becomes much larger than k^. The principal effect is a reduced influence of 
kz on frequency or applied field and thus the longitudinal modes form a continuous band. When d/L=0.1, 
the difference between the first and the fifth longitudinal mode is of O.SkOe at fixed frequency (20GHz). 
The same difference is only of O.OSkOe when d/L=0.02, while the difference in the radial modes is much 
larger. 

2.1 Spatial distribution of edge modes 

The terms which contain the derivatives of the demagnetizing field that appear in the Eq.® have an 
increased influence near the edges of the cylinder where their variation is most important. These terms 
will provide a correction to the frequencies of the edge modes. To quantify this correction we numerically 
calculated the frequencies of the edge modes for a Nickel cyhnder of aspect ratio d/L=0.1 (L=300nm) 
with and without the additional terms at fixed magnetic field. For example, without correction we obtain 



dm 

^r— ±p X m = 



(10) 



min{r,z) ~ Jn{qr)cos{k^z) 



(11) 
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Figure 2: (a) Spatial distribution of an edge mode calculated from the analytical expression for a cylinder 
with diameter-to-length ratio of 0.1 (L=300nm) in a false color code. Red (positive) represents a permitted 
value for the wavenumber. (b) Edge mode frequencies determined from the analytical expression (circles) 
and from micromagnetic simulation (rectangles) at different DC magnetic fields. 

a frequency of 13.79GHz for the edge mode at l.SkOe, a 48% difference from the 9.3GHz value obtained 
with the correction terms. The frequencies obtained for the edge modes by the analytical formula with 
the correction terms are represented in the Figl2](b) along with the values obtained from micromagnetic 
simulations discussed in the next section. These values are similar which imply that the correction terms 
are necessary to calculate exactly the frequencies of edge modes in finite size nanowires. 

In Figl2](a), we show the simulated 2D spatial distribution of the edge mode near surface for a DC 
magnetic field of 2.9kOe. The spatial distribution is obtained from the permitted values of the wavenumber 
near surface (positive). Our method is the following: we first numerically determine the frequencies of the 
edge modes near boundaries assuming that it does not propagate into the center of the cylinder. Next, we 
fix the frequency at the value obtained in the first step (at a given DC magnetic field) and we calculate 
the spatial distribution of the allowed values for the wavenumber in the edge region. From the plot we 
observe that the edge mode is not restricted to the actual surface but goes underneath the surface in 
parabolic fashion. 

3 Micromagnetic simulation 

In order to understand the mode structure and to compare the analytical results of the previous section, a 
micromagnetic simulation was conducted with the nmag package |T6] for a cylinder with an aspect ration 
d/L=0.1 (L=300nm). The cylinder was discretized with a cell size of 5nm. First, the magnetization was 
relaxed to equilibrium using a large damping parameter q:=0.5. The equilibrium configuration was then 
excited adiabatically with a small 20GHz rf magnetic field in the perpendicular direction, while a constant 
magnetic field was applied along the cylinder axis. The values used for the parameters are the same as 
used in the analytical simulation and are those of Ni: D=2x IQ-^Gcm^, 7=188.5GHz/T (g factor of 2.15) 
and M5=480 emu/cm^ [8. 17j and a damping parameter of 0.015. We subsequently computed the spatial 
distribution of the oscillation of the dynamic magnetization doing Fourier transforms on a number of cells 
along the z and x axes of the cylinder. For comparison with previous results, we kept fixed the frequency 
of the rf magnetic field and we varied the amplitude of the DC magnetic field. 

To determine the spatial distribution of the longitudinal modes in an cylindrical nanowire we calculated 
the amplitude of oscillation of the magnetization at different magnetic fields. The values of the static 
magnetic field were chosen to correspond to longitudinal modes as calculated with Eq.([8j). In Figl3](a), we 
show the average of the magnetization profiles along the z-axis for one DC magnetic field which should 
correspond to the 5th longitudinal spin wave mode determined from the analytical model. These profiles 
are calculated doing Fourier transforms on the magnetization for each nm in the z direction and in the 
X direction and extracting the values at 20GHz. Basically, we compute the amplitude of the dynamical 
magnetization variation in a xz-plane passing through the center of the cylinder (Figl3](b)) and then we 
average between the curves in the x direction. This corresponds with the average of a contour plot in 
a xz-plane. As observed, the calculated profile corresponds well with the one expected at this applied 
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Figure 3: (a) Average magnetization profile along z axis of a cylinder with diameter-to-length ratio of 0.1 
(L=300nm). The average was calculated from 30 magnetization profiles taken each nm in the x direction 
in the plane y=0; (b) 2D color projection in the xz plane of the 30 magnetization profiles. The maximal 
amplitude of the dynamical magnetization (red) precess at 20GHz with a DC magnetic field applied 
(2.9kOe) corresponding to the 5th longitudinal mode. 




Figure 4: 2D color projection in the xz plane of two edge modes of a cylinder with diameter-to-length 
ratio of 0.1 (L=300nm) at two different DC magnetic fields: (a) 1.5kOe and (b) 2.9kOe. 

external field from the analytical model. The spins which form the longitudinal mode precess are mainly 
in the effective length Az. Another mode appear at the edges of the cylinder where the spins precess at 
a lower frequency. 

In FiglH the 2D spatial distribution of the magnetization is shown corresponding to the edge modes 
obtained at two different DC magnetic fields. The xz plane shown passes through the center of the cylinder, 
rotating it 360° around its axis will provide the 3D profile of the mode. The panel (b) corresponds to an 
edge mode with spins precessing at a lower frequency of 13.25GHz (DC field 2.9kOe) than the longitudinal 
mode shown in the Figl3](b) where the spins precess at 20GHz. The spatial distribution of edge modes 
comes into agreement with that calculated in the previous section (Figl2]^a)). The frequencies of the edge 
modes obtained from the microwave simulation are shown in the Figl2]^b). The values are similar with 
those calculated with the Eq.(|8]). 

In FigEl a snapshot of the 3D profile of the y component of magnetization is shown after 30ps. We 
clearly observe the variation of the dynamical magnetization along the z direction (long axis). The mode 
is 3D longitudinal mode as expected even though some hybridization effects may exists, the standing 
waves picture being simplistic jl^. The edge mode is clearly seen, it corresponds to the spins rotating at 
the surface (in plane) and precessing at lower frequency than the bulk modes. The edge mode is localized 
and it is always present due to the magnetic field inhomogeneity at the boundaries. The longitudinal 
modes which precess at 20GIIz decay in time, becoming zero after 1.5ns. 

4 Discussion and conclusion 

The spatial distribution of longitudinal and edge modes in cylinders of moderate aspect ratio was cal- 
culated in this paper. It provides new insight on the spatial distribution of spin waves in cylindrical 
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Figure 5: 3D profile of the evolution of the y component of magnetization is shown at a given instant in 
time (30ps). The applied DC magnetic field is 2.9kOe and the frequency of the rf field is 20GHz. 

nanowires. We used a 2D analytical model and obtained the spin wave modes frequencies and profiles for 
the edge modes. The validity of the 2D model was tested through comparison with 3D micromagnetic 
results. The values of the resonances agree quantitatively although the profiles of the edge modes from 
3D micromagnetics are not exactly parabolic and seem to penetrate on a longer distance beneath the 
surface (30-50nm). The difference between the two results are probably due to the assumptions and 
approximations that where made: in the analytical model the magnetization is considered uniform in the 
whole cylinder and a correction in the z direction is obtained but the variation of the radial demagnetizing 
field is not considered (which can be an useful extension of the theory) while in the 3D micromagnetics 
we break the cylinder into small tetrahedral cells (discretization length of 4.8nm) and consider that the 
magnetization is continuous and vary linearly and then the demagnetizing field is computed. 

The nanowires used in the simulation have finite diameter-to-length ratio (d/L=0.1) and thus the 
longitudinal modes are well separated and their spatial distribution was resolved. We tried to resolve the 
spatial distribution of longitudinal modes for a cylinder with d/L=0.02 without success, as the modes are 
hybridized (closer in frequency) and form a quasi-continuum band. Furthermore, as the demagnetizing 
field is inhomogeneous it confines the spin waves at boundaries of the cylinder: the edge modes are always 
present in finite aspect ratio cylinders. These modes can be observed in experiments. Until now, only 
modes of an ensemble of nanowire have been investigated [T7] . To measure the spectrum of an individual 
nanowires a particularly well suited technique is the (Ferro) Magnetic Resonance Force Microscopy [I9l [20] . 
This local probe technique can validate our results. 
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Appendix 

Noting m+ = rux-l-imy et m_ = mx — irriy (similarly to the transverse right and left handed magnetization 
density operators) we can rewrite Eq.Q as: 



4^ - Nr 



, d . d . , d . d . 



ON,, 
dz 



(12) 



Applying the operators (^ — z^)m+ and ( 



A. 

dx 



to Eq.®, it can be easly shown that: 
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+ [uj' - j\H, - DV^f]Ms^ + (4^ - Nrrh^M, 

x{H,-DV^)\/l^m = (13) 

with = V2-^. 

Replacing the scalar potential into Eq. (fT3|l we obtain Eq.([8]). 
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